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Abstract — A  novel  and  efficient  method  for  calibrating  a  sensor 
array  with  position  uncertainties  is  proposed  in  this  paper.  The 
method  is  based  on  two  non-disjoint  sources  in  unknown 
directions  and  three  carry-on  instrumental  sensors.  It  can  be 
applied  to  arbitrary  array  geometries  including  linear  arrays. 
Besides,  no  small  position  error  assumption  is  made,  which  is 
always  an  essential  prerequisite  for  many  existing  array  shape 
calibration  techniques.  The  new  method  achieves  a  favorable 
array  shape  calibration  just  using  a  one-dimensional  search, 
with  no  high-dimensional  nonlinear  search  and  convergence 
burden  involved.  It  is  also  possible  to  extend  the  proposed  idea 
to  tackle  the  problem  of  direction  dependent  gain  and  phase 
uncertainties.  Simulation  results  are  provided  to  demonstrate 
the  effectiveness  and  behavior  of  the  proposed  method. 

I.  Introduction 

Since  their  introduction,  high-resolution  direction  finding 
algorithms  such  as  MUSIC  and  ML  have  received  significant 
attention.  This  is  due  to  their  potential  ability  to  resolve 
sources  separated  by  less  than  one  standard  beamwidth  of  the 
receiving  array,  unlike  the  conventional  Fourier-based 
direction  finding  procedures.  Yet  despite  this  potential 
advantage  offered  by  high-resolution  methods,  their 
application  to  real  systems  has  been  very  limited.  One  of  the 
main  reasons  for  this  situation  is  the  practical  difficulty 
associated  with  calibrating  the  array  manifold,  since  the 
performance  of  these  methods  depends  strongly  on  the 
accuracy  of  the  array  manifold.  However,  in  practice,  the 
actual  array  manifold  always  differs  from  the  nominal  array 
manifold  due  to  sensor  position  errors,  gain  and  phase 
perturbations,  mutual  coupling,  etc.  The  presence  of  these 
unknown  calibration  errors  is  the  major  factor  limiting  the 
performance  of  the  high-resolution  methods  in  practical 
direction  finding  system.  Hence  to  achieve  high-resolution 
performance,  array  calibration  is  always  necessary.  The 
primary  interest  of  this  paper  is  focused  on  the  sensor 
position  errors,  although  the  method  proposed  here  can  also 
be  easily  extended  to  the  case  of  gain  and  phase  uncertainties. 

Sensor  position  errors  induce  the  direction  dependent 
sensor  phase  errors  to  the  array  manifold  so  as  to  exert 
detrimental  effects  on  the  performance  of  the  high-resolution 
direction  finding  algorithms.  Various  techniques  have  been 
developed  in  the  literature  to  circumvent  the  problem  of 
sensor  position  errors.  They  are  either  active  and  need 
calibrating  sources  in  known  directions  or  passive  and  rely 
upon  the  sources  present  in  the  field  to  achieve 
self-calibration.  To  date,  however  most  of  these  techniques 


do  not  work  perfectly  in  the  sense  that  they  are  unable  to 
always  acquire  satisfactory  array  shape  calibration  due  to 
convergence  burden  of  multimodal  nonlinear  search  [l]-[5], 
or  if  they  can,  they  are  either  too  costly  to  implement  due  to 
the  need  of  auxiliary  calibrating  sources  in  known  directions 
[6]-[8]  or  some  certain  pathological  array-source  geometries 
disable  them  at  all  [9]-[  10],  Besides  almost  all  these  above 
array  shape  calibrating  techniques  assume  that  the  position 
perturbations  are  relatively  small  deviations  from  the 
nominal  positions  and  thus  a  first  order  approximation  to  the 
perturbed  array  response  vector  is  often  used  to  simplify  the 
estimation  procedures.  However,  many  simulation  results  in 
[4]  show  that  these  techniques  fail  under  even  moderate 
perturbation  errors.  The  motivation  of  this  paper  is  to  attempt 
to  suggest  an  efficient  and  relatively  inexpensive  and 
practical  scheme  for  the  array  shape  calibration. 

The  scheme  proposed  here  needs  two  non-disjoint  sources 
in  unknown  directions.  It  relaxes  the  small  error  assumption 
and  search  convergence  burden.  The  ambiguity  problem  of 
linear  array  identified  in  [9]-[  1 0]  can  also  be  mitigated.  The 
only  price  paid  for  above  merits  is  that  three  carry-on 
instrumental  sensors  are  needed  to  work  as  coordinate 
reference  and  at  the  same  time  introduce  some  more  degrees 
of  freedom  to  tackle  the  identifiability  problem  associated 
with  the  linear  array. 

II.  Problem  Formulation 

Consider  an  array  of  K  sensors  of  arbitrary  geometry 
impinged  by  M=  2  uncorrelated  sources  from  far  field  in 

unknown  directions  at  8  =  \81,9~>Y  ■  The  signal  waveforms 
are  assumed  to  be  narrowband  of  known  center  frequency. 
The  actual  positions  of  these  K  sensors  differ  from  their 
nominal  positions.  In  addition,  we  have  three  carry-on 
instrumental  sensors,  whose  positions  are  assumed  precisely 
known  and  one  of  which  is  chosen  as  the  origin  of  the 
coordinates.  As  a  result,  an  array  of/V=K+3  sensors  is  formed. 
The  complex  envelope  of  the  noise-corrupted  N  sensors  array 
output  vector  X [t)  may  be  written  as  (1): 

X(t)=  A(0)S{t)+  N{t)  t  =  1,2,  •••Z  (1) 

where  A(?)  is  a  Mx  1  signal  vector,  N(t)  is  a  Nx  1  noise 
vector.  It  is  assumed  that  the  signals  and  noises  are  stationary, 
zero  mean  uncorrelated  Gaussian  random  processes  and 
further,  the  noises  are  both  spatially  and  temporally  white 
with  variance  a1  .Array  manifold  matrix  A[0 )  is  NXM 
matrix  whose  columns  are  the  steering  vectors.  In  the 
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presence  of  sensor  position  uncertainties,  the  N*1  steering 
vector  W(()m )  can  be  modeled  as  (2): 

w(0rn)=r(0m)a(0m)  m  =  1,2  (2) 

where  a(dm )  is  the  ideal  steering  vector  corresponding  to 
array  nominal  positions  and  can  be  expressed  as  (3)-(4)  : 

a(em)=  1,  exp 

*mn  =[xn>yn][sin{Om)cos{emY  n  =  l,2,  —  N  (4) 
1  is  the  wavelength  of  the  signal,  \xn  yn\  is  the  nominal 
co-ordinates  of  the  nth  sensor  with  respect  to  the  reference 
sensor.  Further,  r{()m)  is  a  N*N  complex  diagonal  matrix 
whose  nnth  entries  are  the  angularly  dependent  phase 
distortion  induced  by  the  nth  sensor  position  error  and  it  can 
be  written  as  (5)-(6): 

r(em)=  \,exp^-j^-Arm^,-,exp^-j^-AxmN^  m  =  1,2  (5) 

4Tmn=[dxn,4y„][sin(0m)cos(0m)Y  n  =  l,2,---N  (6) 
where  [Axn  Ayn  ]are  the  position  disturbance  associated  with 
the  nth  sensor. 

The  array  covariance  matrix  and  its  eigendecomposition 
are  expressed  as  follows: 

R  =  E[x{t)XH(t\  =ARsAh  +cj2I  (7) 

M  N  tQ\ 

«  =  S  V/«f  +  IVa"  =  esases  +  enanen  (») 

z=l  i=M+ 1 

where  {l,- ;  i  =  1,2,  •  •  • ,  N;  Xt  >  XM  }  and  { e ; ;  i  =  1,2,  •  •  • ,  N}  are 
ordered  eigenvalues  and  corresponding  eigenvectors  of  R 
respectively.  The  signal  subspace  and  noise  subspace  of R  are 
respectively  the  ranges  of  the  matrices: 

Es=[ei,e2,---eM]  (9) 

En  =  \?  M+l’e  M+2’"' ’e  n\  (10) 

The  problem  of  interest  here  is  as  follows:  given  L  array 
snapshots  x{t)  t  =  1,2,  •  •  • ,  L  ,  estimate  the  unknown  DOAs  of 
sources,  as  well  as  the  unknown  sensor  position  uncertainties 
of  the  K  sensors,  \_Axn  Ayn  ]  n  =  4,5,-  •  • N  (N  =  K  +  3) . 

III.  Algorithm  Description 

Since  we  have  3  instrumental  sensors  with  no  sensor 
position  uncertainties,  the  first  3  diagonal  entries  of  matrix 
r(d )  are  all  reduced  to  1.  The  ideal  steer  vector  a(d)  and 
phase  distortion  matrix  r(d )  can  be  partitioned  as  follows: 

a(9)=[al(9)  4(4  (11) 

r(0)  =  diag[lu3  [vecdfT  2)]'\  (12) 

where  the  3  x  1  vector  ax  ($)  and  K  x  1  vector  «2  ($)  are  formed 
from  the  part  of  the  elements  of  a{0)  corresponding  to  the 
instrumental  sensors  and  the  position-disturbed  sensors 
respectively.  Similarly,  the  diagonal  entries  of  KxK  diagonal 
matrix  r2  {&)  consist  of  the  unknown  direction  dependent 
phase  uncertainties  induced  by  the  corresponding 

position-disturbed  sensors.  Then  the  steer  vector  tV(d)  can 
be  reformulated  as  follows: 


m  =  1,2  ^ 


w(o)=r(e)a(p) 

ax{6)  0ixK  1 
_0Kx  1  diag[a2(o)]\yecd{r2{o))\ 

=  a{e)s{d) 

where  a{6)  is  a  tfxjf  +  l)  matrix  while  S[d)  is  a 
(Ai  +  l)xl  vector.  t//ag[v]  denotes  a  diagonal  matrix  whose 
diagonal  entries  is  formed  from  the  elements  of  vector  v  and 
vecd[A ]  denotes  a  column  vector  where  the  diagonal 

elements  of  A  form  the  vector. 

The  underlying  basis  for  subspace -based  DOA  estimation 
algorithms  is  the  orthogonality  between  the  noise  subspace 
and  signal  subspace  of  array  covariance  matrix  R.  which 
means  that 

WH{e)ENE“w{d)=  0  (14) 

Substitution  of  (13)  into  (14)  yields  (15-17): 

dH{d)iH{d)ENE^a{d)8{d)^  0  (15) 

sH  (e)Q{e)s{e)= o  (16) 

Q{d)  =  aH{d)ENE^a{d)  (17) 

where  Q[d)  is  a  [K  +  l)x(*  + 1)  Hermitian  matrix.  Since  the 
d(d)^0  ,  (16)  means  that  the  matrix  Q[6 )  is  singular, 
i.e.,  ran  k[Q{e)]  <K  + 1 .  Note  that  under  the  condition  that 
array  manifold  \w(0) :  -  k/2  <6  <  tt/2}  is  unambiguous,  the 
matrix  Q(O)  is  singular  or  rank  reduction  if  and  only  if 
the  9  =  Oj  i  =  1,2  ,  since  the  dimension  of  signal  subspace  off? 
is  2.  Based  on  this  idea,  we  develop  a  DOA  estimator  as  (18) 
or  (19)  and  a  array  shape  calibration  algorithm  as  (20)-(26): 

f)  =  argmax - r-/  A  (18) 

e  1  0(9) 


9  =  arg  max  - 


det  Q9 


s{9)=e,mn[Q[9)\  with  emin(l)  =  l  (20) 

Q(9)  =  aH(9)ENE^a(9)  (21) 

vecd  (f2(^|=[j(2)43),- ,S{K  +  l)J  (22) 

AX  =  [Ax4  Ax5---Axn]T  (24) 

AY  =  [Ay4Ay5...AyN]T  (25) 

p(0m )  =  ~^-angle[vecd(f2 i(Jm )]  m  =  1,2  (26) 

2k 

where  EN  denotes  the  finite  sample  estimate  of  noise 
subspace  .  Tmin  [(g(6>)]  is  the  smallest  eigenvalue  of  the 
matrix  Q{o) ,  emm  [gfe)]  is  the  eigenvector  corresponding  to 


the  smallest  eigenvalue  of  the  matrix  Q\9)  and  the 

det  [(?($)]  is  the  determinant  of  the  matrix  Q(9) . 

From  (18)-(26),  we  observe  that  the  DOA  estimation  and 
array  shape  calibration  can  be  achieved  simultaneously  just 


637 


using  an  one-dimensional  search  over  FOV(Field  of  view)  of 
array.  Besides  with  uniform  linear  array,  due  to  the 
Vandermonde  structure  of  the  ideal  steering  vector,  a 
polynomial  rooting  with  degree  of  4  can  also  be  utilized  in 
light  of  the  idea  behinds  ROOT-MUSIC  algorithm. 

Although  the  small  position  error  assumption  is  not  made 
in  the  above  shape  calibration  algorithm,  it  is  always  assumed 
that  the  position  uncertainties  are  not  large  to  such  an  extent 
that  the  phase  delay  ambiguity  identified  in  [11]  is  present. 

Finally,  it  is  obviously  that  the  idea  proposed  here  can 
easily  be  extended  to  array  calibration  in  the  presence  of 
direction  dependent  gain  and  phase  calibration. 

IV.  Simulation  Results 

In  this  section,  simulation  results  are  presented  to  illustrate 
the  performance  of  the  new  algorithm.  Simulations  are 
carried  out  for  a  nominal  uniformly  linear  array  of  1 6  sensors 
with  one  half-wavelength  inter-sensor  spacing.  Three 
instrumental  sensors  are  added  to  form  a  nominal  uniformly 
linear  array  of  19  sensors.  The  corresponding  reference 
coordinate  system,  nominal  array  geometry  and  perturbed 
array  geometry  are  shown  in  the  fig.  1.  The  actual  sensor 
positions  are  arbitrarily  fixed  and  allowed  to  randomly  vary 
from  the  nominal  sensor  positions  within  the  range  ±  0.5  A  in 
Y-coordinates  and  ±  0.25  A  in  X-coordinates.  Two 
narrowband  uncorrelated  sources  with  equal  power  impinge 
on  the  array,  from  the  far  filed,  at  distinct 
directions  30°  and  40°  w.r.t  the  broadside  of  array.  The 
SNR=20dB  is  defined  as  the  ratio  of  each  signal  power  to  the 
noise  power  at  each  sensor.  200  snapshots  are  used  to 
estimate  the  array  covariance  matrices.  The  number  of 
sources  is  assumed  known. 

Fig. 2  shows  the  spatial  spectra  obtained  form  the  new 
algorithm  (18).  In  table  1-2,  we  also  demonstrate  the 
estimated  value  and  real  value  of  the  sensor  positions. 

From  the  results  presented  above  and  many  other 
simulations  with  similarly  favorable  results,  it  can  be 
concluded  that  the  array  shape  calibration  algorithm 
proposed  here  provides  us  an  efficient  and  relatively 
inexpensive  array  shape  calibration  scheme. 


[2]  AJ.  Weiss  and  B.  Friedlander,  “Array  shape  calibration  using 
eigenstructure  methods,”  Signal  Processing,  vol.  22,  no. 3,  pp.  251-258, 
Elsevier  Science  Publishers,  1991. 

[3]  J.  H.  Hong,  “Genetic  approach  to  bearing  estimation  with  sensor 
location  uncertainties”.  IEE  Electronic  letter  Vol. 29  no.23, 
pp.2013-2014.  Nov  1993. 

[4]  B.  P.  Flanagan  and  K.  L.  Bell,  “Improved  array  self-calibration  with 
large  sensor  position  errors  for  closed  space  sources”.  Proc.  2000 
Sensor  Array  and  Multichannel  workshop,  Cambridge,  MA,  March 
2000,  pp.484-488. 

[5]  C.  R.  Wan,  J.  T.  Goh  and  C.  S.  Chia,” Array  calibration  based  on 
Inverse  subspace  fitting”,  IEEE  2001  Conference  and  Exhibition  on 
OCEANS  Vol  4,  pp.2468-2491. 

[6]  B.P.Ng,  S.Wee  “Array  shape  calibration  using  source  in  known 
locations’Woc  ICCS/ISITA  ’92  Sigapore  pp. 836-840. 

[7]  B.  C.  Ng  A.  Nehorai  “Optimum  active  array  shape  calibation”  Proc. 
25th  Asilomar  Conference  on  signal, system  and  computers  1991 
pp. 893-897. 

[8]  Ming  Zhang  and  Zhao-Da  Zhu  “Array  shape  calibation  using  sources 
in  known  directions”.  Aerospace  and  electronics  Conference  NAECON 
1993  pp. 70-73. 

[9]  Y.  Rockah  and  P.  M.  Schultheiss,  “Array  shape  calibration  using 
sources  in  unknown  locations-Part  I:  Far  field  sources”,  IEEE  Trans 
ASSP  Vol.35,  no.3  pp286-317  Mar  1987. 

[10]  Y.  Rockah  and  P.  M.  Schultheiss,  “Array  shape  calibration  using 
sources  in  unknown  locations-Part  II:  Near  field  sources  and  estimator 
implementation”,  IEEE  Trans  ASSP  Vol.35,  no.6,  pp724-735,  Mar 
1987. 

[11]  J.T.-H.Lo  and  S.  L.  Marple,  “Observability  conditions  for  multiple 
signal  directon  finding  and  array  sensor  localization”.  IEEE  trans  SP 
Vol.40,  no.  1 1 ,  pp.264 1  -265,  Nov.  1 992 


nominal  array  shape  and  perturbed  array  shape 


perturbed  array  shape 
•  nominal  array  shape 


L  instrument  sensers 


o 


o 


3  4  5  6 

X(wavelength) 


Fig.  1.  Reference  array  coordinate  system 


V.  Conclusion 

In  this  paper,  we  propose  an  efficient  array  shape 
calibration  algorithm  by  using  two  non-disjoint  sources  in 
unknown  direction  and  three  instrumental  sensors.  It  can  be 
applied  to  arbitrary  array  geometries  including  linear  arrays. 
The  new  method  is  computationally  abstractive  and  relatively 
inexpensive.  Besides,  small  position  error  assumption  is 
relaxed  to  meet  the  need  of  array  shape  calibration  in  the 
presence  of  large  position  errors.  Without  any  modification, 
the  idea  proposed  here  can  also  be  extended  to  the  array 
calibration  in  the  presence  of  direction  dependent  gain  and 
phase  uncertainties. 

References 

[1]  A.  J.  Weiss  and  B.  Friedlander,  “Array  shape  calibration  using  sources 
in  unknown  locations-A  maximum  likelihood  approach”,  IEEE  Trans 
ASSP  Vol. 37  no.  12,  pp  1958- 1966,  Dec.  1989. 


Fig.  2.  Spatial  spectrum  acquired  with  new  method 


638 


TABLE  1 

Estimated  X-Coordinates 


Coordinate 

s 

Nominal 

Actual 

Estimated 

X4 

1.5 

1.3454 

1.3354 

X5 

2.0 

2.1719 

2.1783 

X6 

2.5 

2.3370 

2.3403 

X7 

3.0 

2.8354 

2.8450 

X8 

3.5 

3.7471 

3.7465 

X9 

4.0 

3.9699 

3.9840 

X10 

4.5 

4.4200 

4.4321 

Xll 

5.0 

4.9071 

4.9029 

X12 

5.5 

5.4325 

5.4262 

X13 

6.0 

5.9466 

5.9390 

X14 

6.5 

6.5458 

6.5479 

X15 

7.0 

6.8099 

6.8093 

X16 

7.5 

7.2691 

7.2664 

X17 

8.0 

7.9793 

7.9755 

X18 

8.5 

8.6849 

8.6864 

X19 

9.0 

9.2171 

9.2150 

TABLE  2 

ESTIMATED  Y-COORDINATES 


Coordinate 

s 

Nominal 

Actual 

Estimated 

Y4 

0 

-0.2356 

-0.2257 

Y5 

0 

-0.3397 

-0.3426 

Y6 

0 

0.3729 

0.3721 

Y7 

0 

-0.2621 

-0.2687 

Y8 

0 

0.1458 

0.1447 

Y9 

0 

0.4669 

0.4560 

Y10 

0 

0.1649 

0.1547 

Yll 

0 

0.3704 

0.3711 

Y12 

0 

-0.4901 

-0.4884 

Y13 

0 

-0.3630 

-0.3602 

Y14 

0 

0.3188 

0.3151 

Y15 

0 

-0.0698 

-0.0706 

Y16 

0 

0.3903 

0.3921 

Y17 

0 

0.2349 

0.2389 

Y18 

0 

0.1873 

0.1888 

Y19 

0 

-0.1539 

-0.1521 
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